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The noisy Burgers equation in one spatial dimension is analyzed by means of the Martin- Siggia- 
Rose technique in functional form. In a canonical formulation the morphology and scaling behavior 
are accessed by mean of a principle of least action in the asymptotic non-perturbative weak noise 
limit. The ensuing coupled saddle point field equations for the local slope and noise fields, replacing 
the noisy Burgers equation, are solved yielding nonlinear localized soliton solutions and extended 
linear diffusive mode solutions, describing the morphology of a growing interface. The canonical 
formalism and the principle of least action also associate momentum, energy, and action with a 
soliton-diffusive mode configuration and thus provides a selection criterion for the noise-induced 
fluctuations. In a "quantum mechanical" representation of the path integral the noise fluctuations, 
corresponding to different paths in the path integral, are interpreted as "quantum fluctuations" and 
the growth morphology represented by a Landau-type quasi-particle gas of "quantum solitons" with 
gapless dispersion E oc P^/^ and "quantum diffusive modes" with a gap in the spectrum. Finally, 
the scaling properties are dicussed from a heuristic point of view in terms of a "quantum spectral 
representation" for the slope correlations. The dynamic eponent z = 3/2 is given by the gapless 
soliton dispersion law, whereas the roughness exponent ^ = 1/2 follows from a regularity property 
of the form factor in the spectral representation. A heuristic expression for the scaling function 
is given by spectral representation and has a form similar to the probability distribution for Levy 
flights with index z. 
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I. INTRODUCTION 

This is the second of a series of papers where we analyze the Burgers equation in one spatial dimension with the 
purpose of modelling the growth of an interface; for a brief account we refer to In the first paper, denoted in the 
following by A , we investigated the noiseless Burgers equation [^|-^ in terms of its nonlinear soliton or shock wave 
excitations and linear diffusive modes. In the present paper we address our main objective, namely the noisy Burgers 
equation in one spatial dimension [0. This equation has the form 

— =v\/'^u + \u\/u + \/ri , (1.1) 
at 

where v is a damping constant or viscosity and A a nonlinear coupling strength. The equation is driven by a conserved 
white noise term, Vr/, where 77 has a Gaussian distribution and is short-range correlated in space according to 

(r/(a;, i)ry(a;', t')) = M{x - x')S{t - t') . (1.2) 

In the context of modelling a growing interface the Kardar-Parisi-Zhang equation (KPZ) for the height field h [|j 

^^^V^h+^mr+rj, (1.3) 

is equivalent to the Burgers equation by means of the relationships 

u^Vh (1.4) 
'udx, (1.5) 
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that is the Burgers equation governs the dynamics of the local slope of the interface. In Fig. 1 we have sketched the 
growth morphology in terms of the height and slope fields. 

The substantial conce ptu al problems encountered in nonequilibrium physics are in many ways embodied in the 
Burgers-KPZ equations (|l . l| ) and (|l.3|) which describe the self-affine growth of an interface subject to annealed noise 
arising from fluctuations in the drive or in the environment P ITsl. Interestingly, the Burgers-KPZ equations are 
also encountered in a variety of other problems such as randomly stirred fluids [[7|, dissipative transport in a driven 
lattice gas [p^-[l^ jthe propagation of flame fronts , the sinc-Gordon equation l22| , and magnetic flux lines in 

superconductors |2^. Furthermore, by means of the Cole- Hopf transformation ||2^ , |25|| the Burgers-KPZ equations 
are also related to the problem of a directed polymer |^,^ or a quantum particle in a random medium ]2^ , |29| and 
thus to the theory of spin glasses |]30|~^ . 

In contrast to the case of the noiseless Burgers equation discussed in A where the slope field eventually relaxes due 
to the dissipative term z^V^m, unless energy is supplied to the system at the boundaries, the noisy Burgers equation 
( |1 . 1[ ) describes an open nonlinear dissipative system driven into a stationary state with random energy input at a 
short wave length scale provided by the conserved noise, V77. In the stationary regime the equation thus describes 
time-independent stochastic self-affine roughened growth. In the linear case for A = the Burgers equation reduces 
to the noise-driven Edwards- Wilkinson equation (EW) ||3|] 

^yV^u + Vrj- (1.6) 

at 

here for the slope field u. Owing to the absence of the nonlinear growth term AmVu the cascade in wave number space 
is absent and the correlations, probability distributions, and scaling properties are easy to derive ]34| . Furthermore, 
unlike the Burgers case the EW equation, being compatible with a fluctuation-dissipation theorem, actually describes 
the dynamic fluctuations in an equilibrium state with temperature A/2i^ (in units such that fc^ — 1) and as a 
consequence does not provide a proper description of a growing interface. On the other hand, the presence of the 
nonlinear growth term \uVu in Eq. (^]^) renders it much more complicated and much richer. The term filters the 
input noise V77 and gives rise to interactions between different wave number components leading to a cascade which 
changes both the scaling properties and the probability distributions from the linear EW case. 

The Burgers-KPZ equations owing to their simple form accompanied by their very complex behavior have served 
as paradigms in the theory of driven and disordered systems and have been studied intensively p|-p^. One set of 
issues which have been addressed are the scaling properties [^-^. According to the dynamic scaling hypothesis 
|35lp8[ supported by numerical simulations and the fixed point structure of a renormalization group scaling analysis, 
I I^QTthe slope field u is statistically scale invariant in the sense that the self-affine rescaled u'(a;, t) — b^'^'^^^^u{bx, b^t) 
is statistically equivalent to u{x,t), where 6 is a scale parameter. More precisely, the scaling hypothesis implies the 
following dynamical scaling form for the slope correlation function in the stationary regime [|7|,p|j3q-|4Q| : 



{uix, t)u{x: t')) ^\x- («-i)/(|t - t'\/\x - xT ) . (1.7) 

The scaling behavior in the long wave length-low frequency limit is thus governed by two scaling dimensions: i) the 
roughness or wandering exponent ^, characterizing the slope correlations for a stationary profile and ii) the dynamic 
exponent z, describing the temporal scaling in the stationary regime |^l[]. The slope field u has the scaling dimension 
1 — C- For large w the scaling function /(w) oc u'~^^^~'>^/^; for small w f{x) oc const. 

Two properties determine the scaling exponents, namely a scaling law and an effective fluctuation-dissipation 
theorem. Like the noiseless or deterministic Burgers equation discussed in A, the noisy equation is also invariant 
under a nonlinear Galilean transformation ||^,|^ 

x—i- X — Xuot (1-8) 
u + uq . (1-9) 

Since the nonlinear coupling strength A here enters as a structural constant of the Galilean symmetry group it 
transforms trivially under a scaling transformation and we infer the scaling law p|-pi][| 



C + z = 2, (1.10) 



relating C and z. Furthermore, noting that the stationary Fokker-Planck equation for the Burgers equation (LI) is 
solved by a Gaussian distribution |0,O,p^ 



P{u) cx exp 



'a 



dxu^ 



(1.11) 
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i nde pendent of A it follows that u is an independent random variable and that the height variable h according to Eq. 
(1.5) perfor ms ra ndom walk, corresponding to the roughness exponent C = 1/2, also in the linear EW case. From the 
scaling law (1.10) we subsequently obtain the dynamic exponent z = 3/2. In the linear EW case C = 1/2 and z — 2, 
characteristic of diffusion. In Table 1 we have summarized the exponents for the EW and Burgers/KPZ universality 
classes. 



Universality class 


c 


z 


EW 


1/2 


2 


Burgers-KPZ 


1/2 


3/2 



Table 1. Exponents and universality classes 

The standard tool used in the analysis of the scaling properties of nonlinear Langevin equations of the type in 
Eqs. (l.l-OI) is the dynamic renormalization group (DRG) method |^. The method is based on an expansion 
in powers of the nonlinear couplings, an expansion in the noise strength, arising from the noise contractions when 
implementing the statistical average, combined with an infinitesimal momentum shell integration in the short wave 
length limit, corresponding to the scaling transformation. The method operates in variable dimension d and typically 
identifies critical dimensions separating regions where infrared convergent perturbation theory holds yielding mean 
field behavior from regions with infrared divergent expansions. Here the DRG allows for an organisation of the 
divergent terms and yields renormalization group equations for the effective parameters in the theory in terms of an 
epsilon expansion about the critical dimension . In powers of epsilon the DRG thus yields expression for the critical 
exponents and information about the scaling functions. 

For the Burgers-KPZ equations the expansion is in powers of the nonlinear coupling strength A and the noise 
strength A. The critical dimension is dc = 2, also following from simple power counting. Below dc — 2 there appear 
two renormalization group fixed points: An unstable Gaussian fixed point, corresponding to vanishing coupling 
strength, describing a smooth interface governed by the EW equation, and a stable strong coupling fixed point, 
characterizing a rough interface. The exponents assume non-trivial val ues fo r all A 0. In c? = 1 an effective 
fluctuation-dissipation theorem equivalent to the Gaussian form in Eq. (1.11) is operative, and together with the 
Galilean invariance, implying trivial scaling of A, the renormalization group equations yield the exponents in Table 

1 associated with the infrared stable non-trivial strong coupling fixed point [7[PJ4^-^. Unlike the case of static 
and dynamic critical phenomena, where the renormalization group methods have proven very successful ]35|j4^ , the 
situation in the case of nonequilibrium phenomena exemplified here by the Burgers-KPZ equations has proven more 
difficult and despite extensive efforts based on DRG calculations to higher loop order and mode coupling 
approaches p5[ | the physics of the strong coupling fixed point in rf = 1 still remains elusive. 

In a recent letter, denoted in the following by L [ ^Tf , we approached the strong coupling fixed point behavior from 
the point of view of the mapping of the Burgers equation onto an equivalent solid-on-solid or driven lattice gas model 
|^,|8|, which furthermore maps onto a discrete spin 1/2 chain model [ ^9|j5C| ]. The quantum spin chain approach has 
been proposed in |]ll| , ^ , and considered further in , on the basis of the equivalence between the Liouville operator 
in the Master equation describing the evolution of the one-dimensional driven lattice gas, or the equivalent lattice 
interface solid-on-solid growth model, and a non-Hermitian spin 1/2 Hamiltonian. The quantum chain model has 
been treated by means of Bethe-Ansatz methods JlT| , ^l| -^ and the dynamic exponent z = 3/2 obtained from the 
finite size mass gap scaling. In L we pushed the analysis further and constructing a harmonic oscillator representation 
valid for large spin in combination with a continuum limit we derived a Hamiltonian description and a set of coupled 
field equations of motion for the spin field, corresponding to the slope u, and a conjugate "azimuthal" angle field 
replacing the noise. The field equations admit spin wave solutions, corresponding to the linear diffusive modes and, 
more importantly, nonlinear localized soliton solutions, describing the growing steps in the original KPZ equation or 
the solitons or shocks in the Burgers equation. We also derived the soliton dispersion law and after a quasi-classical 
quantization identified in a heuristic manner the elementary excitations of the theory. From the dispersion law we 
deduced the dynamic exponent z = 3/2, characteristic of the zero temperature fixed point of the "quantum theory". 
The picture that emerged from our analysis was that of a dilute quasi-particle gas of nonlinear soliton modes yielding 

2 = 3/2 and a superposed spin wave gas, corresponding to z = 2, the dynamic exponent for the linear case. In L 
we also briefly discussed the operator algebra associated with the Hamiltonian representation and derived the fleld 
equations by means of a canonical representation of the Fokker-Planck equation for the equivalent Burgers equation. 
Whereas the Bethe-Ansatz investigations by their nature are restricted to special values of the coupling strength, 
corresponding to the fully asymmetric exclusion model ji^ , our analysis is valid for general coupling strength and 
thus constitutes an extension of the Bethe-Ansatz method to the general case of a continuum field theory. The 
analysis in L was in many respects incomplete and preliminary but it did indicate that the strong coupling flxed point 
behavior is intrinsically associated with the soliton modes in the Burgers equation since they both provide aspects of 
the growth morphology and also, independently, yield the dynamic exponent. 
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Here we present a unified approach to the noisy Burgers based on the Martin-Siggia-Rose technique in functional 
form |p6H6l|. This method supersedes the analysis in L and does not make use of the mapping to a spin chain via 
a solid-on-solid model and therefore the implicit assumption of persistent universality classes under the mappings. 
Clearly, the different formulations are basically equivalent arising as they do from the same basic stochastic growth 
problem. The equivalence also indirectly demonstrates that the scaling properties of the different models, i.e., the 
solid-on-solid model, the spin chain, and the field formulation, fall in the same KPZ-Burgers universality class. The 
present functional path integral method which also can be seen as a generalization of the stationary distribution (|l.ll|) 




to the time-dependent nonlinear case provides a many-body description of the morphology of a growing interface in 
terms of soliton excitations and does also give insight into the scaling behavior. Below we highlight some of our 
results. 

• The path integral formulation yields a compact description of the noisy Burgers equation and provides expres- 
sions for the probability distributions and correlation functions. Reformulated as a canonical Feynman-type 
phase space path integral the approach allows for a principle of least action. Hence the weight of the different 
paths or interface configurations, corresponding to the noise-induced interface fluctuations, contributing to the 
path integral are controlled by an effective action. The role of the effective Planck constant is here played by 
the noise correlation strength A. The action in the path integral thus plays the same role for the dynamical 
configurations as the Hamiltonian in the Boltzmann factor for the static configurations in equilibrium statistical 
mechanics. 

• In the asymptotic weak noise limit the principle of least action implies that the dominant configurations arising 
from the solutions of the saddle point field equations correspond to a a dilute nonlinear soliton gas with super- 
posed linear diffusive modes. The canonical formulation and the principle of least action furthermore allow a 
dynamical description and associate energy, momentum, and action with the solitons and the diffusive modes. 

• The path integral formulation permits a "quantum mechanical" interpretation in terms of an underlying non- 
Hermitian relaxational "quantum mechanics" or "quantum field theory" . The noise-induced fluctuations here 
correspond to "quantum fluctuations" and the fluctuating growth morphology is described by a Landau-type 
quasi-particle gas of nonlinear "quantum solitons" and linear "quantum diffusive modes" . In the height field this 
corresponds to a morphology of growing steps with superposed linear modes. The "quantum soliton" dispersion 
law is gapless and characterized by an exponent 3/2; the "quantum diffusive mode" dispersion law is quadratic 
with a gap in the spectrum proportional to the soliton amplitude. 

• In the present formulation the scaling properties associated with the "zero temperature" fixed point in the 
underlying "quantum field theory" follow as a by-product from the soliton and diffusive mode dispersion laws and 
the spectral representation of the correlations. The dominant excitation in the long wave length-low frequency 
limit identifies the relevant universality class. The present many-body formulation yields the known exponents. 
The dynamic exponents z = 3/2 and z = 2 are associated with the soliton and diffusive mode dispersion laws, 
respectively, whereas the roughness exponent C = 1/2 follows from a regularity property of the form factor in 
the spectral representation. The many-body formulation also explains the robustness of the roughness exponent 
under a change of universality class and provides a heuristic expression for the scaling function which has the 
same structure as the probability distribution for Levy flights. 

• From a field theoretical point of view we identify the noise strength A as the effective small parameter. Further- 
more, the fundamental probability distribution or path integral has an essential singularity for A = 0. Hence 
our approach is based on a non-perturbative saddle point or steepest descent approximation to the path integral. 
We believe that it is in this respect that the dynamic renormalization group method based on an expansion in 
A and in the noise contraction A fails to access the strong coupling fixed point. 

The path integral representation of the noisy Burgers equation presented here is equivalent to a full-fledged one 
dimensional non-Hermitian non-Lagrangian field theory and requires for its detailed analysis some advanced field 
theoretical techniques and methods from quantum chaos. In the present context we choose, however, a somewhat 
heuristic approach to the path integral in order to elucidate the emerging simple physical picture of a growing interface. 
This approach then also serves as a tutorial introduction to the field theoretical treatment to be presented elsewhere. 

The present paper is organized in the following way. In section II we discuss the simple case of the linear Edwards- 
Wilkinson equation, mainly in order to emphasize the non-perturbative nature of the noise as regards the stationary 
driven regime. Since the soliton modes in the noisy Burgers equation turn out to be of crucial importance in under- 
standing the morphology and scaling properties, we summarize in section HI the results obtained in A concerning the 
solitons and diffusive modes in the noiseless Burgers equation. In section IV we set up the path integral formulation 
for the noisy Burgers equation in terms of the Martin-Siggia-Rose techniques in functional form. In section V we 
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perform a shift transformation of the path integral to a canonical Feynman path integral form and discuss the canon- 
ical structure and the associated symmetry algebra. Section VI is devoted to an asymptotic weak noise saddle point 
approximation and to the derivation of the deterministic coupled field equations replacing the Burgers equation. In 
section VII we solve the field equations and derive nonlinear soliton and linear diffusive mode solutions. In section 
VIII we discuss the dynamics of the solitons following from the principle of least action. The dominating morphology 
of a stochastically growing interface can be interpreted in terms of a dilute soliton gas; this aspect is discussed in 
some detail in section IX. The fluctuation spectrum about the soliton solutions is basically given by the path integral, 
however, in section X we take a heuristic point of view and discuss the fluctuations as "quantum fluctuations" in 
the underlying non-Hermitian "quantum field theory" . Section XI is devoted to a discussion of the scaling properties 
and universality classes on the basis of the "elementary excitations" in the "quantum description" . We also present 
a heuristic expression for the scaling function. Finally in section XII we present a discussion and a conclusion. 



II. THE LINEAR EDWARDS - WILKINSON EQUATION - THE ROLE OF NOISE 

Here we review the properties of the linear case described by the noise-driven Edwards- Wilkinson equation 
in particular in order to elucidate the role of the noise. For the slope field this equation is given by 

^^^^^u + Vt] (2.1) 
ot 



with the noise r/ correlated according to Eq. (|L 

{tj{x, t)r]{x',t')) = AS{x - x')6{t - t') . (2.2) 
The equation ( |2.lD has the form of a conservation law 

l^-V, (.3, 

with current 

j = -vVu - 77 . (2.4) 

We note that with average vanishing Vu at the boundaries the conservation law implies that the average off-set in 
the height, J Whdx, is conserved. 



In wave number space, u(k,t) = J dccexp (—ikx)u{x,t), and solving Eq. (2.1 as an initial value problem averaging 
over the noise according to Eq. (pj) we obtain for the slope correlations 

{u{k,t)u{~k,t')) = [(u(/c,0)u(-A:,0))i- A/2]exp[-(t + t')i^fc2] + A/2exp[-|i-tVA:2] _ (^2.b) 

Here {■ ■ ■)i denotes an average over initial values which is assumed independent of the noise average (•••). The basic 
time scale is set by the wave number dependent lifetime T(fc) = l/vk"^] which diverges in the long wave length limit 
/c 0, characteristic of a conserved hydrodynamical mode. We note that at short times compared to T(k), which 
sets the time scale for the transient regime, {u{k, t)u{—k, t')) is non-stationary and depends on the initial correlations, 
whereas at long times t, t' 3> T(k) the correlations enter a stationary, time reversal invariant regime and depends only 
on |< — For vanishing initial slope, u(fc,0) = 0, we obtain in particular the mean square slope fluctuations 

(l"(A^,0l'>-y[l-exp[-2t/r(fc)]] (2.6) 
which approaches the saturation value A/ 2 for t ^ rik) 



More precisely, it follows from Eq. (2.5) that for fixed t — t' the transient term can be neglected at times greater 



than a characteristic crossover time tco of order 

U l/(i^fc') log (1/A) , (2.7) 

depending also on the noise strength A. This time thus defines the onset of the stationary regime. For t ^ T{k), tco 
noise-induced fluctuations built up and the mean square slope fluctuations approach the constant value A/2. For 
A 0, tco ~^ 00 and the system never enters the stationary regime. 
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The "elementary excitation" is the diffusive mode u{k,t) oc exp(±i/fc^t). In frequency space u{k,uj) = 
J dtexp {iujt)u{k,t) and the slope correlation function assumes the Lorentzian diffusive form, characteristic of a 
hydrodynamical mode, 



{u{k,uj)u{~k,—Lo)) = 



Ak^ 



(2.8) 



with diffusive poles at = zLii'k'^, a strength given by A/i^ and a line width vk^. We note that in the stationary 
regime both the decaying and growing modes, u oc exp (ii^/c^t), contribute to the stationary correlations. Time 
reversal invariance is thus induced from the microscopic reversibility of the noise-driven system. In the transient 
regime for t <C T{k) the initial conditions enter and we must choose the solution propagating forward in time, 
u oc exp {—vkHV in order to satisfy causality. 

From Eq. (2.8) we also obtain the scaling function 



fH 



(A/2i^)(47riy)~i/2u;-i/2 exp [-l/Ai^w] 



(2.9) 



in accordance with the general form in Eq. ( |l.7| ) yielding the EW exponents in Table 1 defining the EW universality 
class. For large w f(w) ^ w'^^"^; for small w f(w) — > but with an essential singularity for w = 0. In frequency-wave 
number space the scaling form is 



(u(fc, cj)u(— fc, — ti>)) = fc '^g[ijj/k'^) 



and we directly infer the scaling function 



9{w) = 



A 



(2.10) 



(2.11) 



In Fig. 2 we have shown the slope correlation function and the scaling functions / and g in the EW case. 

In contrast to the noisy Burgers equation, the EW equation does not provide a proper description of a growing 
interface. This is seen by expressing Eq. (2.1) in the form 



du 
'dt 



,5F_ 
6u 



where the effective free energy is given by 



F = - / dxu^ 



(2.12) 



(2.13) 



Using the fluctuation-dissipation theorem to relate the noise strength to an effective temperature T it then follows 
that the EW equation describes time-dependent fluctuations in an equili brium system with temperature T — A/2iy 
and with an equilibrium distribution given by the Boltzmann factor Eq. (1.11), i.e., 



P{u) (X exp 



(2.14) 



We already here in the linear case note that the noise strength A seems to play a special role. Whereas A enters 
linearly in the correlation function (uu)(fc , m), the limit of vanishing noise strength, A — > 0, appears as an essential 
singularity in the stationary distribution (2.14). Since the distribution P{u), appropriately generalized to the time- 
dependent case, is the generator for the correlation function (uu) and higher commulants, it is clearly the fundamental 
object and the role of the noise strength A as a non-perturbative parameter an important observation. More precisely 
the point is the following: Whereas the damping constant v together with the relevant wave number k defines the time 
scale for the transient regime where the system has memory and evolves forward in time in an irreversible manner, 
the presence of the noise is essential in order for the system to leave the transient regime at all and to enter the 
stationary regime where t he s ystem is time reversal invariant, as for example reflected in the evenness in lu in the 
slope correlation function (2.8). In the absence of the noise the system simply decays owing to dissipation unless it is 



driven by deterministic currents at the boundaries. Imposing the noise and driving the system stochastically is thus 
a sm^wZar process, as reflected mathematically by the essential singularity in the distribution ( 2.14 ). 

Although the above observation of the non-perturbative role of the noise strength A is a trivial statement in the 
linear case where is just reflects the structure of the Boltzmann factor, we will later show that in a more complete 
theory of a growing interface, described by the nonlinear noisy Burgers equation, it is essential to take into account 
non-perturbative contributions in the noise strength A. 
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III. THE SOLITON MODE IN THE NOISELESS BURGERS EQUATION 



It turns out that the sohton excitation in the noiseless Burgers equation when properly generahzed play an important 
role in the understanding of the growth morphology and strong coupling behavior of the noisy Burgers equation. In A 
we discussed in some detail the soliton and diffusive mode solutions in the noiseless Burgers equation and performed 
a linear stability analysis. Here we briefly summarize those aspects of the analysis in A which will be of importance 
in the discussion of the noisy case. 

The noiseless or deterministic Burgers equation has the form [BH6I 



— = v\7'^u + \u\7u. (3.1) 

and is a nonlinear diffusive evolution equation with a linear term controlled by the damping or viscosity v and a 
nonlinear mode coupling term characterized by A. In the context of fluid motion the nonlinear term gives rise to 
convection as in the Navier Stokes equation; for an interface the term corresponds to a slope dependent growth. 

Under time reversal t —t —t and the transformation u — )■ —u the equation is invariant provided v — > —v. This 
indicates that the linear diffusive term and the nonlinear convective or growth term play completely different roles. 
The diffusive term is intrinsically irreversible whereas the growth term corresponding to a mode coupling leads to a 
cascade in wave number space and a generates genuine transient growth. The transformation t ^ —t \s absorbed in 
u — *■ —u or, alternatively, A — > —A, corresponding to a change of growth direction. We also note that the equation 
is invariant under the parity transformation x —x provided u —u. This feature is related to the presence of 
a single spatial derivative in the growth term and implies that the equation only supports solitons or shocks with 
one parity, that is right hand solitons. Finally, the Burgers equation is invariant under the Galilean symmetry group 
( pr8[]L9| ), that is a Galilean boost to a frame moving with velocity Auq is absorbed by a shift in the slope field by uq. 



The irreversible and diffusive structure of Eq. (3.1) implies that an initial disturbance eventually decays owing to 
the damping term vV^u. In the linear case the slope field decays by simple diffusion u{x, t) esc exp {—vkF't) exp (±.ikx) 
as discussed in section II. In the presence of the nonlinear mode coupling term the equation also supports localized 
soliton or kink profiles ||6^-|65|] with given parity. In the static case the symmetric positive parity or right hand soliton 
has the form 

u{x) = w+ tanh [ks{x — xa)] (3-2) 
ks = Xu+/2v . (3.3) 

We have introduced the characteristic wave number k^ setting the inverse length scale associated with the static soliton, 
xq denotes the center of mass position. The width of the soliton is of order l/k^ and depends on the amplitude W-f. 
In the inviscid limit u or for strong coupling A — > cx), the wave number fc^ — > oo and the soliton reduces to a 
sharp shock wave discontinuity. 



Boosting the static soliton in Eq. (3.2) to a finite propagation velocity v and at the same time shifting u we obtain. 



denoting the right and left boundary values by u±, the soliton solution 



u[x, t) = ^ 1 ^ tanh [— («+ — u-){x — vt ~ xo)\ (3-4) 



with velocity v given by the soliton condition 



2v . . 

u+ + ti_ = -— . (3.5) 
A 



We note that the soliton condition (3.5) is consistent with the fundamental nonlinear Galilean invariance and remains 
invariant under the transformation: v ^ v + Amq and u± — s- u± — uq. Also, unlike the case for the Lorentz invariant 
(j)'^ and sine-Gordon evolution equations [6^ ], the propagation velocity in the present case is tied to the amplitude 
boundary values of the soliton - a feature of the nonlinear Galilean invariance. In Fig. 3 we have depicted the right 



hand soliton solution given by Eq. (3.4) and the associated height field h 



In the linear case for A = case the diffusive modes with dispersion 

Lol ^ -ivk^ (3.6) 

"exhaust" the spectrum of relaxational modes. For A 7^ the soliton profile acts as a reflectionless Bargman potential 
giving rise to a bound state at zero frequency, corresponding to the translation mode of the soliton - the Goldstone 
mode restoring the broken translational invariance, and a band of phase-shifted diffusive scattering modes [p5|. The 
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resulting change of density of states is in accordance with Levinson's theorem in that the potential traps a bound 
state and depletes the continuum of one state. In the presence of the soliton the diffusive modes furthermore develop 
a gap in the spectrum of uik as depicted in Fig. 4, 

ujk = -ii^ik^ + kl) . (3.7) 

An asymptotic analysis of the noiseless Burgers equation in the inviscid limit v [ p6[ shows that an initial 
configuration breaks up into a "gas" of propagating and coalescing kinks connected by ramp solutions of the form u cx 
const — x/Xt. This allows for the following qualitative picture of the transient time evolution: Although the nonlinear 
mode coupling term is incompatible with a proper superposition principle we can still along the lines of the evolution 
of integrable one dimensional evolution equations envisage that an initial configuration "contains" a number of 
right hand solitons connected by ramps. In the course of time the solitons propagate and coalesce. Superposed on 
the soliton gas is a gas of phase-shifted diffusive modes. As discussed in A the gap in the diffusive spectrum can be 
associated with the current flowing towards the center of the solitons. The damping of the configuration predominantly 
takes place at the center of the soliton where u varies rapidly thus enhancing the damping term i/V^u. We also note 
that only parity breaking right hand solitons are generated in the noiseless Burgers equation. In Fig. 5 we have shown 
the transient evolution of the slope field and the associated height field. 



IV. PATH INTEGRAL REPRESENTATION OF THE NOISY BURGERS EQUATION 



In this section we begin the analysis of the noisy Burgers equation. In our discussion in section II of the linear EW 
equation we noticed that the noise strength A enters in a non-perturbative way in the stationary distribution in Eq. 
( |DI1 ). Whereas this, of course, is a trivial observation in the linear case since the EW equation describes fluctuations 
in equilibrium and A oc T, that is the singularity structure is the same as the low T limit of the Boltzmann factor 
exp (—E/T), the presence of the nonlinear mode coupling growth term in the noisy Burgers equation renders the 
situation much more subtle. We are now dealing with an intrinsically nonequilibrium situation. The noise drives the 
system into a far-from-equilibrium stationary state and equilibrium statistical mechanics does not apply. On the other 
hand, from our study of the noiseless Burgers equation, we have learned that the soliton excitations play an important 
role in the dynamics of the morphology of a growing interface and is a direct signature of the nonlinearity. The issue 
facing us is then how to include both the non-perturbtive aspects of the noise and the nonlinear soliton structure 
in a consistent way. It turns out that the functional formulation of the Martin-Siggia-Rose techniques provides the 
appropriate formal and practical language for such an approach ||5^-|6l[ . 

Our starting point is the noisy Burgers equation (1.1) for the fluctuating slope field u, i.e.. 



du 

at 



which has the structure of conserved nonlinear Langevin equation with current 

A 2 

J — — vu u — 77 . 



(4.1) 



(4.2) 



For the noise we assume a Gaussian distribution 



P(?7) oc exp 



1 

'2A 



dxdtri{xty 



(4.3) 



where 77 is correlated according to Eq. (1.2), i.e., 

{7]{X, t)7]{x', t')) = M{x - x')5{t - t') 



(4.4) 



Unlike the transient relaxation of an initial value configuration described by the deterministic Burgers equation, the 
noisy Burgers equati on is driven continuously by the conserved noise V77, corresponding to a fluctuating component of 
the current j in Eq. ( |4.2| ) . Energy is fed into the system via the noise and dissipated by the linear damping term. The 
nonlinear mode coupling gives rise to a cascade in wave number space corresponding to "dissipative structures" in the 
growth morphology. This mechanism changes the probability distributions and associated correlatio ns (m oments). 



scaling exponents, and scaling functions from the EW case in section II. In other words, the equation (4.1 
nonlinear box which transforms the input noise V77 to an output slope field u. 



acts as a 
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In the Martin-Siggia-Rose techniques the probability distribution for the slope field P(u) and the correlations (uu) 
are conveniently derived from an effective partition function or generator 67 



Z{n) = ( exp 



i / dxdtu{x,t)^{x,t) 



(4.5) 



Here fi{x, t) is a generalized chemical potential or external conjugate field coupling to the slope field u{x, t) and (• • 
denotes an average over the input noise r/, implementing the nonlinear stochastic relationship provided by the Burgers 
equation ( |4.l| ). In terms of Z we have for example the probability distribution P{u) = {5{u — u(x, t)))^, 



P{u{x,t)) 
and the correlation function 

{u{x,t)u{x',t')) 



exp 



—i / dxdtu{x , t) fj,{x , t) 



Z{^{x,t)) 



5Z{ii) 



Sfi{x, t)Sfi{x' , t') 



(4.6) 



(4.7) 



higher moments are derived in a likewise manner. In order to incorporate the nonlinear constraint imposed by the 
Burgers equation we insert the identity 



(4.8) 



in the partition function Z{fi)] for a first order evolution equation one can show that causality implies that the 
Jacobian relating du to du/dt equals unity Finally, e xpo nentiating the delta function constraint in Eq. ( |4.8| ) 

and averaging over the noise distribution according to Eq. (4.3) we obtain 



Z{^) — j JJ^ dwdpexp [iG] exp i j dxdtu^L 



where the effective functional G is given by 

G = / dxdt 



p(^ ~ ^^^"^ ^ AuVii) + ^A(Vp)^ 



(4.9) 



(4.10) 



Th e path or functional integral (4.9) with G given by Eq. (4.1C) effectively replaces the stochastic Burgers equation 
(4.1). The path integral is deterministic and the noise rj is replaced by the different configurations or paths contributing 
to Z . In this sense Z is an effective partition function for the dynamical problem and G an effective Hamiltonian, 
analogous to the Hamiltonian in the partition function ^ = '^'^P i^H/T) in equilibrium statistical mechanics. We 
also note that the transcription of the Burgers equation to a path integral leads to t he a ppearance of an additional 
noise field p, arising from the exponentiation of the delta function constraint in Eq. (4.8) |5^-^, and replacing the 
stochastic noise in Eq. ( [4.1[ ). 

Since the path integral formulation provides a field theoretical framework allowing for functional and diagrammatic 
techniques, Feynman rules, skeleton graphs. Ward identities, etc., it is mostly used in order to generate perturbation 
expansions in powers of the the nonlinear coupling \uS/u p^- ^ , |67| , |69| -|7l]] . It is, however, worthwhile noting that 
such a fiel d th eoretic expansion has precisely the same structure as the one produced by directly iterating the Burgers 
equation (4.1) in powers of AmVu and averaging over the noise term by term according to Eq. (|4.4|). 



Also, corroborating our remarks in section II, we notice from the structure of the path integral ( [4.9| - [4.10D that the 
noise strength A appears as a singular parameter in the sense that A — s- gives rise to the singular delta function 
constraint for the Burgers equation. This limit is, however, much more transparent when we express Z in a "canonical 
form" . 



CANONICAL TRANSFORMATION TO A HAMILTONIAN FORM 

INTEGRAL - SYMMETRIES 



PHASE SPACE PATH 



By inspection of the path i ntegra l in Eqs. (4.9-4.10) we notice that it has the same structure as the usual phase 
space Feynman path integral [|67|]7^j7^ as regards the kinetic term pdu/dt in F but that otherwise p and u do not 



9 



appear in a canonical combination. This situation can, however, be remedied by performing a simple complex shift 
of the noise variable p 



(5.1) 



in Eqs. (OfLlO). Assuming that the path integral operates in a space time LT box, i.e., |a;| < L/2 and \t\ < T/2, and 
imposing periodic or vanishing boundary conditions for u and (p in order to eliminate total derivatives, the partition 
function Z^jj,) can be expressed as 



Z{fi) = const J Y\_dud(pexp ^'^'^ ^^P * J 



i / dxdtufi 



where the action S is given by the canonical form 

S — ( dxdt 

with the complex Hamiltonian density 



- n{u,ip) 



2 ^ 



(5.2) 



(5.3) 



(5.4) 



The Hamiltonian density consists of two terms: A relaxational or irreversible harmonic component, — «(z^/2)[(Vu)^ + 
(V<^)^], corresponding to the diffusive aspects of a growing interface, i.e., the linear damping, and a nonlinear reversible 
mode coupling component, (A/2)u^V(p, associated with the drive A. 

One feature of the transcription of the noisy Burgers equation to a canonical path integral form is that the effective 
Hamiltonian density (5.4) driving the dynamics of the system is in general complex. This particular aspect was also 
encountered in the treatment in L where the growth term in the spin chain Hamiltonian turned out to be complex. 
We also notice that the doubling of variables, i.e., the replacement of the stochastic noise rj by an additional noise 
field ^p in the path integral, was also encountered in the treatment in L in the canonical oscillator representation of 
the spin variables. 

It is here instructive to compare the above path integral for the relaxational growth dynamics of the Burgers 
equation with the usual phase space path integral formulation in quantum mechanics or quantum field theory | |72| , [73t . 
Here the partition function has the form 



Z = 



xt 



with the classical action 



S 



dxdt 



dq 



nip,q) 



(5.5) 



(5.6) 



where p and q are c onsidere d canonical ly conju gate variables and 7i(p, q) the usual classical Hamiltonian density. 

Comparing Eqs. (5.5-5^) with Eqs. (5.2-^T^) it is evident that the structures of the two path integral formulations 
are quite similar and we are led to identify the noise strength A/v with an "effective" Planck constant and the 
Hamiltonian density Ti. as the generator of the dynamics. The classical limit thus corresponds to the weak noise limit 
A and in analogy with the quasi-classical or WKB approximation in quantum mechanics, A — > 0, constitutes 
a singular limit in accordance with our previous remarks. The partition function Z with the action S given by Eqs. 
( ^.2[j5.4D t hus constitutes the required generahzation of the stationary distribution P{u) oc exp[— (j//A) J dxu"^] in 
Eq. (1.11) to the time-dependent case |75|. By com parison we furthermore conclude that the slope field u and the 
noise field ip, replacing the Gaussian noise in Eq. ( |4.l|) , are canonically conjugate momentum and coordinate variables 
satisfying the Poisson bracket algebra iTq] 



{u{x), fix')} — d{x — x') 



and that the Hamiltonian or energy 

H= f dxH 



dx 



[(V^.)2 + (V^)2]+ Vv^ 



(5.7) 



(5., 
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is the generator of time translations according to the equations of motion 



du ^ ^ , , 

— ={il,u} (5.9) 

^^{H,^}. (5.10) 

Drawing on the mechanical analogue the momentum P, the generator of translations in space, is also easily identified 
from the basic transformation properties, 

Vu = {P,u} (5.11) 

V^ = {P,(/^}, (5.12) 

and it follows that 



P = J dxg (5.13) 
g = mV(^ , (5.14) 

where g is the momentum density. 

In order to elucidate the canonical structure of the path integral ( ^.2[-^.4D and the analogy with the usual phase space 
Feynman path integral we have generated a complex Hamiltonian ( |5.4| ). Note, however, that by formally rotating the 
noise field in phase space if ^ iip the Hamiltonian and the action become purely imaginary leading to a real path 
integral. 

The symmetries discussed in the context of the c^uantum spin chain representation in L are also easily recovered 
here. Noting that H is invariant under a constant shift of the noise field, ip ^ (p + (po, we infer that the integrated 
slope field 



M J dxu (5.15) 

i.e., the total off-set of the height field, h = J dxu, across the interface, is a constant of motion, 

{iJ,M} = 0. (5.16) 

This is consistent with the local conservation law following from the structure of the Burgers equation, but is here a 
consequence of the structure of the path integral. The invariance under a shift of ip is equivalent to the invariance of 
the Burgers equation under a shift of the noise rj in the noise term V?/. Similarly, under a constant shift of the slope 
field u u + uq, we have, introducing the momentum density g, 7i — s- 7i + Xuog + (A/2)iiQV.g, or since the last term 
is a total derivative, H ^ H + AwqP, corresponding to an associated Galilean transformation with velocity — Awq. 
For the integrated noise field 

$ = J dxip (5.17) 

we thus obtain the Poisson bracket algebra 

{iJ,$} = AP (5.18) 

which together with 

{H,M} ^ {H,P} ^0 (5.19) 
{P, $} = {P, M} = (5.20) 

and 

{$,Af} = P (5.21) 



defines the symmetry algebra. We note again that the nonlinear coupling strength A enters the Poisson bracket ( 5.1§| ) 
and thus is a structural constant of the symmetry group. 
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We finally wish to comment on the properties of the path integral in Eqs. (5.2-5.4) under time reversal t — > —t. By 
construction the path integral applies at late times compared to any initial time to, defining the initial value of the 
slope configuration uq- This implies that the noise in the Burgers equation has driven the system into a stationary 
time regime and that the transients associated with uq have died out. In the linear case for A = 0, we note by 
inspection that the path integral is invariant under the combined operation i — > — i and ip — > — implying that 
the slope correlations are not only stationary but also invariant under time reversal. This is in agreement with the 
analysis of the noisy EW equation in section II where we obtained {uu){k,uj) = Ak'^/[uj'^ + {i/k'^)'^], implying that 
{uu){x,t) depends on \t\. This is consistent with the description of an equilibrium interface and is just an expression 
of microscopic reversibility. In the presence of the drive for A 7^ the path integral is invariant under the com bined 
transformation t —> —t, Lp — and A — > —A showing that the term (A/2)w^V</3 in the Hamiltonian (5.8) gives 
rise to a proper growth direction thereby breaking time reversal invariance, that is we are dealing with a genuine 
nonequilibrium phenomena. 



VI. FIELD EQUATIONS IN THE WEAK NOISE LIMIT - SADDLE POINT APPROXIMATION 



The basic structure of the path integral (5.2-5.4) is illustrated by the simple one dimensional integral, 



/(A) = / duexp 



i-^S{u) 



exp [ifiu\ , 



(6.1) 



where A is the small parameter (the noise strength). In the limit A — > the integral /(A) is approximated by a steepest 



descent calculation which amounts to an expansion of S{u) about an extremum uq, S{u) 
and a subsequent calculation of a Gaussian integral. For small A we then obtain 



/(A) = exp [i—S{uo)] exp [i/xuo] exp 



2S"{uo) 



-27riA 



1/2 



S{uo) + ^S"{uo){u - uo) 



(6.2) 



The leading contribution to /(A) is given by exp [zS'(uo)/A] and is thus determined by the extremal value of the 
action. This part, however, goes along with a multiplicative factor, [— 27riA/S"'(uo)]^/^, arising from the Gaussian 
integral sampling the fluctuations about the stationary points; this term is the first in an asymptotic expansion in 
powers of A^/^. We notice that there is an essential singularity for A = 0, signalling the non-perturbative aspects 
of a steepest descent calculation; the result cannot be obtained as a perturbation expansion in powers of A. The 
analysis of the path integral now essentially follows the same procedure but is rendered much more difficult owing to 
the field theoretical phase space structure of the problem. In Fig. 6 we have depicted the principle of a saddle point 
or steepest descent calculation of /(A). 

In the weak noise limit A — > the asymptotically leading contribution to the path integral thus arises from 
configurations or paths (u, if) corresponding to an extremum or stationary point of the action S. Invoking the 
variational condition ^iS* = with respect to independent variations of the slope field u and the canonically conjugate 
noise field ip, Su and Sip, with vanishing variations at the boundaries of the space time LT box, we readily infer, using 
Eq. (5.8), the classical equations of motion 



du SH 
ot du 



Implementing the functional derivation or, equivalently, using the Poisson bracket relation (5.7), we obtain 



du . 2 ^ 

— — —IV\ ip + AUVU 

ot 

-^= +ih'V^u + XuVip 
ot 



(6.3) 
(6.4) 

(6.5) 
(6.6) 



The above coupled field equations (6.5- p.q ) are a fundamental result of the present analysis. They provide a 
deterministic description of the noisy Burgers equation in the asymptotic non-perturbative weak noise limit. The 
equations have the same form as the ones derived in L based on the quasi-classical limit of the quantum spin chain 
representation. Furthermore, the parameter identification is in accordance with the "quantum representation" of the 
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Fokker-Planck equation. As regards the considerations in L this demonstrates that the precise identification of the 
quasi-classical hmit is in fact a weak noise limit in the exact path integral representation of the Burgers equation. 

First of all we observe that the field equation for the slope field u has the form of a conservation law, du/dt = — Vu, 
with current j — —{X/2)u^ + ivW^p. The fluctuating component in the current in the noisy Burgers equation, 
j = i(A/2)u^ — v'S/u — rj, is thus replaced by the noise field ip and admissible solutions must yield an imaginary noise 
field in order to render a real current, corroborating our remarks in the previous section. The field equation for the 
noise field is parametrically coupled to the slope field and in the presence of the coupling A driven by the momentum 
density g = uVip. 

Secondly, we confirm that the field equations are invariant under the nonlinear Galilean transformation (1.8-1.9|), 



and under an arbitrary shift in tp 



u{x,t)^ u{x — XuQt,t) — uq (6-7) 
ip{x,t)^ ip{x — Xuot,t) (6-8) 



(p{x,t) ^ (p{x,t) - ipo . (6.9) 



in accordance with the general discussion of the symmetry algebra and consistent with the symmetry properties of 
the noiseless and noisy Burgers equations. 

One final comment on the classical zero noise limit. We maintain that in the asymptotic non-perturbative weak 
noise limit the coupled field equations provide the correct description of the leading behavior of the noisy Burgers 
equation. In order to obtain the noiseless Burgers equation discussed in section III we must confine the noise field 
strictly to the line ip = iu in {u, ip) phase space in which case both field equations reduce to the noiseless Burgers 



equation (3.1). Note, however, that setting (p = —iu we obtain the noiseless Burgers equation with v replaced by 
—V supporting growing linear modes and a left hand nonlinear soliton mode - the missing modes necessary in order 
to describe the correct morphology in the noise-driven stationary regime. In other words, we anticipate that the 
lines p = zkiu define regions for the stationary steepest descent or saddle point solutions of the field equations. The 
vicinity of these lines correspond to the Gaussian fluctuations about the stationary points, that is the linear diffusive 
modes. This picture will in fact be will be borne out when we turn to an analysis of the fleld equations in the next 
section. In Fig. 7 we have shown the extremal paths, corresponding to the saddle point solutions and the nearby 
paths characterizing the fluctuations in (w, p) phase space. 



VII. SOLITON AND DIFFUSIVE MODE SOLUTIONS OF THE FIELD EQUATIONS 



The field equations (3.5-6.6) for u and p constitute a set of nonlinear coupled partial differential equations. The 
general solution is not known. Unlike the noiseless Burgers equation which can be solved by means of the nonlinear 
Cole-Hopf transformation, similar substitutions do not seem to work for the field equations. Presently, it is not known 
whether the field equations belong to the small class of nonlinear evolution equations which can be integrated partly 
or completely by means of the inverse scattering method and related techniques. We are therefore obliged to choose 
a more pedestrian approach and search for special solutions to the equations |62,p3,Ba]. 



A. Stationary states 

We note that the constant slope-constant noise configurations 

u^uo (7.1) 
ip^ po (7.2) 

are trivial saddle point solutions with vanishing energy, momentum, and action. They form an infinitely degenerate set 
and correspond to the zero-energy aligned ferromagnetic spin states discussed in L. The degenerate stationary slope 
configurations are related by a Galilean transformation and we shall in general choose a state with vanishing slope 
corresponding to a horizontal interface. As regards the noise field we are free to choose it equal to zero. In the phase 
space plot in Fig. 7 the background stationary state or the "vacuum" thus corresponds to the origin (u, p) = (0, 0). 
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B. Linear diffusive modes 



In the linear case for A = the Hamihonian H is harmonic in the fields u and f. The coupled field equations are 
linear, 

^= -^i^VV (7.3) 
^= +iiyV\ , (7.4) 

and describe the weak noise limit of the EW equation. Expanding about the stationary state (uqiVo) = (0,0) the 
equations readily admit the solutions 

u{xt)= Y^[u[+\-'^"'+'''^ + u[-\+'^°'-'''^] (7.5) 

k 

with the quadratic diffusive dispersion law 

col = -iiyk^ . (7.7) 

Since u is rea l we have (u^^"*)* = implying that ip is purely imaginary as discussed above. We also note that the 
solution (L5), unlike the solution of the noiseless EW equation, includes both growing and decaying solutions. As 
discussed earlier this feature is consistent with the time reversal invariance in the stationary regime. 



C. Nonlinear soliton modes 



In order to treat the nonlinear aspects of the field equations we employ the same method as in the analysis of 
the noiseless Burgers equation and look for static solutions. Using the Galilean invariance propagating solutions are 
then easily generated by a transformation to a moving frame accompanied by a shift of u. The static case, du/dt = 
dif/dt = 0, corresponds to the action S = ~ J dxdtT-Hu, ip) and the solutions are given by the stationary points of the 
Hamiltonian H. Multiplying the static field equations by W(p and Vit, respectively, we obtain VipV'^ip + VuV^u — 0, 
or by quadrature, imposing the boundary conditions of vanishing slope, Vu = V(/3 = for x —> ±L/2, the slope 
condition 



(V^)2 + {\/uf = . 

Solving the slope condition we obtain 

parametrized by the parity index /i = ±1, which inserted in the static field equations yields 

^u\7^u + AuVu = . 



(7.8) 
(7.9) 
(7.10) 



This equation has the sam e fo rm as the static limit of the noiseless Burgers equation (3.1) with damping /i^, and we 
obtain the static solution (|3.2[) with v replaced by jiv^ i.e., ±i/, 



u{x) = u+ tanh 



2v 



■{x - Xo) 



(7.11) 



The solution (7.11) has the form of a static, localized, symmetric soliton or kink with amplitude 2|m_|.|, center of mass 
xo, width 2j//(A|u+|), approaching ±/i|M+| for x —> ±L/2. In the limit of vanishing damping, v 0, the soliton 
becomes a sharp discontinuity or shock in the slope field. The static soliton connects two degenerate stationary 
states with slopes ±|u+|. However, unlike the sine-Gordon soliton which is characterized by a topological 

quantum number or the tp^ soliton which connects the two degenerate ground states defined by the double-well 

potential, the Burgers soliton has an arbitrary amplitude corresponding to the infinitely degenerate stationary 
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states. Furthermore, we note the interesting fact that unhke the case of the noiseless Burgers equation, where we only 
have a single right hand soliton mode, corresponding to /i = +1, the broken reflection or parity symmetry is restored 
in the noisy case. The noise drives the interface into a stationary state and in the process excites both right and left 
hand solitons. This mechanism in the nonlinear case is equivalent to the excitation of both growing and decaying 
diffusive modes in the EW case. 

The static soliton is a special configuration connecting stationary states with opposite slopes. However, since the 
underlying field equations are invariant under the Galilean symmetry group ( |6.7| - |6.8[ ) it is an easy task to construct a 
propagating soliton solutions by means of a Galilean transformation. Similar to the discussion of the noiseless Burgers 
equation we obtain, introducing the boundary values u± for x — > L/2, the soliton condition (3.5), i.e.. 



and the moving soliton solution 



u{xt) 



■ tanh 



2v 



4v' 



U^\{x — Vt — Xq) 



(7.12) 



(7.13) 



connecting stationary states with slopes m+ and w_; note that we have absorbed the parity index in the sign of 
M_. The soliton has center of mass a;o , propagates with velocity w = — A(m++m_)/2, has the width 4z//(A|m+—m_), 
and amplitude — In the limit of vanishing damping u 0, the inviscid limit, or large drive A, the soliton 
becomes a shock wave, i.e., a discontinuity in the slope field; for finite damping the shock wave front is smoothed 
by dissipation. For vanishing amplitude, u+ — m_ 0, the soliton merges continuously into the stationary state 



Finally, integrating the slope condition (7.8), which by inspection also holds for the propagating soliton, and using 



the invariance property (6.9) in order to set the integration constant equal to zero, we obtain 



If — l/IU 



(7.14) 



yielding the static and propagating soliton solutions for the associated noise field and in accordance with the saddle 
point regions in Fig. 7. In Fig. 8 we have depicted the static solitons and the associated smoothed cusps in the height 
field. 



D. Multi-soliton solutions - Boundary conditions 



In addition to defining the path integral in a finite LT box, we must also specify appropriate boundary conditions 
for the slope and noise fields in accordance with the physical situation. For a growing interface it is convenient to 
assume a horizontal interface at the boundaries equivalent to a vanishing slope field. This corresponds to a vanishing 
deterministic component in the current in Eq. ( 4.2| ) at the boundaries and implies that the interface is only driven 
by the noise. However, since the single soliton solution discussed above connect stationary states with different slopes 
corresponding to a non- vanishing current, we must pair at least two so litons of opposite parity in order to satisfy 
the boundary conditions. By inspection of the field equations (3.5-6.6) we note, however, that a non-overlapping 
two-soliton configuration m*^^) + u^^-* connected by a segment of constant slope is an approximate solution to the 
field equations and therefore corresponds to an extremum of the action in the path integral. The correction term is 
given by A(u^^-' Vu*-^^ + m^^^Vu^^-*) whose contribution to the action we can ignore for non-overlapping well-separated 
solitons. Furthermore, the argument can be generalized to a multi-soliton configuration connected by segments of 
constant slope. It is a well-known feature of path integral instanton or solitons configurations that one in order to 
obtain the correct asymptotic behavior must sum over a gas of non-overlapping instantons or solitons [ f73| , [77[ . The 
situation is the same in the present somewhat more complicated context. In order to satisfy the boundary conditions 
of vanishing slope and to collect all the leading contributions in the asymptotic weak noise limit the structure of the 
path integral implies the formation of a dilute gas of non-overlapping solitons. In Fig. 9 we have shown the case of 
two non-overlapping soliton solutions. 



VIII. DYNAMICS OF SOLITONS - PRINCIPLE OF LEAST ACTION 



It is a fundamental aspect of the canonical form of the path integral for the noisy Burgers equation that it supports 
a principle of least action 1 74 1 . Therefore, unlike the noiseless Burgers equation where there is no underlying canonical 
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structure, the asymptotic weak noise soliton and diffusive mode solutions are derived from a variational principle. 
Furthermore, we can associate an effective action, energy, and momentum with a particular phase space configuration 
or growth pattern. 



The energy density e is generally given by Eq. (5.4). Inserting the slope condition (7 



valid for the soliton 

solutions, the harmonic part of e cancels and the solit on so lutions exclusively contribute to the growth term, i.e.. 



e — {X/2)v?Vip, or in terms of the momentum density ( 5.14 ), e = {\/2)ug. Inserting the soliton constraint (7^) the 
energy density also takes the form 



(A/2)i^u^Vu . 



(8.1) 



We note that the energy density is localized to the position of the soliton where u varies most rapidly. For the soliton 
energy we thus obtain by quadrature in terms of the boundary values m±. 



In a similar manner the soliton momentum (5.14) is given by 

P = 4/1^ [m+ - Ut] 



and finally from Eq. (5.3), using dip/dt = —vVip for the boosted static soliton, the soliton action 

S ^ -T[Pv + E] , 



(8.2) 



(8.3) 



(8.4) 



which also follows from the Galilean invariance of S |74| . 

The purely imagina ry char acter of E, P, and is a feature of our choice of convention in establishing the canonical 
path integral in Eqs. (x2-5^). In order to exploit the formalism of analytical mechanics and the structure of the phase 
space Feynman path integral we have chosen the noise field ip in such a manner that u and ip appear as canonically 
conjugate variables satisfying the usual Poisson bracket (5/7). As discussed earlier this implies that ip has to be purely 
imaginary in the case of the weak noise saddle point solutions. The complex character of E and P is consistent with 
the relaxational and propagating aspects of the soliton modes as also observed in the Fokker-Planck description in L. 
The main properties of E and P in the present dynamical context are that they serve as generators of translations in 
time and space, respectively Jt^ . The nonlinear energy-momentum relationship is characteristic of nonlinear soliton 
solutions [55 1 and is different from the simple E-P relationship encountered in the Lorentz invariant ip'^ or sine-Gordon 
equation [32|. We also note that the damping constant does not enter in the expressions for E and P which only 
depend on the boundary values u± and the drive A. In the weak noise limit the dynamics of soliton solutions is thus 
entirely decoupled from the dynamics of the linear diffusive modes. 

Let us specifically consider a soliton confi gurat ion satisfying the boundary condition of left vanishing slope, i.e., 
U- = for X = —L/2. The soliton condition (7.12) then implies the right boundary value u+ = —2v/X, relating to 
the propagation velocity v. The soliton with positive parity, m+ > 0, propagates left with negative velocity; whereas 
the soliton with opposite parity, < 0, propagates in the forward direction. For E and P wc infer for both parities 



E^'J-^ 
3 A2 

p = -2i- 



S ■ 



A2 



(8.5) 
(8.6) 
(8.7) 



The velocity v = — Aw-|-/2 characterizes the kinematics of the soliton and is related to the amplitude whereas E and 
P determine the transformation properties. Eliminating the velocity we obtain the soliton dispersion law 



Ef 



(8.8) 



where sign(ImP) — — sigm;. We note that the nonlinear localized soliton excitation has a qualitatively different 
dispersion law from the linear extended diffusive mode dispersion law u — —ii/k^. They are both gapless modes 
but the exponents are different. The consequences of this aspect on the spectrum and scaling properties will be 
investigated later when we consider the fluctuations in more detail, but we already note here that the change in the 
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exponent which can be identified with the dynamic exponent z is related to the different universahty classes for the 
EW and Burgers cases. 

Since the energy and momentum densities are localized at the solito n p os itio ns it f ollow s that they are additive 
quantities for a multi-soliton configuration and the general expressions (8.2),(^!^), and (8^) allow us to evaluate the 
total energy, momentum, and action for an arbitrary configuration constructed from well-separated non-overlapping 
solitons. 



IX. A GROWING INTERFACE AS A DILUTE SOLITON GAS 



We are now in position to present a coherent picture of the morphology of a statistically driven growing interface. 
In the weak noise limit A — > the principle of least action which operates in the present context implies that the 
stationary points in (u, ip) phase space correspond to solitons and multi-soliton configurations connected by segments 
of constant slope. In addition there will be superposed diffusive modes. In the nonlinear case the soliton configurations 
determine the dominant features of the growth morphology and will be considered here. The superposed diffusive 
modes will be discussed in the next section. 

The weight of a particular soliton configuration in the path integral is given by the action which is an additive 
quantity for a dilute gas of solitons. The soliton configurations are assumed to be excited with respect to a stationary 
state of vanishing slope, i.e., a horizontal interface, and are furthermore determined by imposing periodic boundary 
conditions at x — ±L/2; we remark that fixed boundary conditions are inconsistent with a soliton configuration 
moving across the system. We also note that unlike the transient properties of the noiseless Burgers equation which 
are described by a gas of right hand solitons connected by ramps, corresponding to a transient height profile composed 
of smoothed cusps connected by convex parabolic segments as shown in Fig. 5 , the stationary state of the noisy 
Burgers equation is characterized by a gas of both right and left hand solitons connected by pieces of constant slope. 
The noise thus radically changes the growth morphology of the Burgers equation. The noise stochastically modifies 
the transient regime by exciting solitons of both parities which thus describe the morphology of the slope field in the 
stationary nonequilibrium state. The situation is similar to the case of the noise-driven damped sine-Gordon equation 
|f78| , [79| where the noise also excites nonlinear soliton modes. 

We now proceed to discuss the morphology of a growing interfa ce in terms of solit ons in t he slope field u. The basic 
"building block" is the static soliton configuration given by Eq. (7.11). From Eqs 
has vanishing momentum P = 0, energy E = i(A/3)|u-|-p, and action S — — i(A/3)T|u 
By integration the height profile, h — J udx, is given by 



2-B.4) it follows that this mode 
independent of its parity. 



h{x) 



2v 



log cosh 



2v 



■{X - Xq) 



(9.1) 



corresponding to a downward and an upward pointing cusp smoothed by the damping constant v. In the limit of 
vanishing ly the cusps become sharp. We also note that the static soliton does not satisfy the boundary conditions of 
vanishing slope. In Fig. 8 we have depicted the slope and height field in the two cases. 

Bo osting the static s oliton in Eq. ( 7.11 ) we obtain a single moving soliton with velocity given by the soliton condition 
( [7.12| ) and the profile ( |7.13| ). In the particular case of a soliton satisfying the left boundary condition u- — 0, we 
obtain from Eq. (7.13) the slope field 



u{xt) 



1 + tanh 



{x ~ Vt — Xq) 



(9.2) 



and by integration the height profile 

h{xt) = 

with propagation velocity 



2i/ 

X |u^ 



log cosh 



Xua 



X\u+ 
Av 



-{x — vt — Xq) 



(9.3) 



(9.4) 



The energy, momentum, and action are given by Eqs. (8.5-8.7). This mode corresponds to the bottom part of an 
ascending step or top part of a descending step in the height field propagating to the left or right, depending on the 
sign of u+. The configurations are shown in Fig. 10. 
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In order to describe a moving step in the height profile we pa ir two well-separated non-overlapping solitons with 
equal amplitude and opposite parity. The soliton condition ( 7.12 ) then implies that they move in the same direction 
with the same velocity. In this case the slope and height fields have the form 



tanh-^— !^(a; 



vt 



h{xt) = 

A \u+\ 
\u+ 



log 



4iy 

cosh {\\u^\/Av){x 



Xi) — tanh 
vt - 



- Xi) 



cosh (A|u+|/4i^)(a; — vt — xi) 



-{x — vt — X2) 



(9.5) 
(9.6) 
(9.7) 



We have here assumed xi <ti X2 for the center of mass coordinates. This configuration corresponds to two co-moving 
solitons moving with velocity v ~ — Au_|_/2 and is equivalent to a moving step in the height profile. The height of 
the step Aft, is given by Aft — uj^(x2 — x\). Imposing periodic boundary conditions for the slope field corresponding 
to a closed ring of length L, this two-soliton mode corresponds to a step in ft moving along the closed ring. At each 
revolution the height field thus increases by Aft and we have a simple gro wth sit uation. For well-separated solitons 
the energy, momentum, and action are additive and we obtain from Eqs. (3.2 -8.4) 



E. 



step 



8 . 



P. 



step 



^ step 



-Ai- 



.v\v\ 



4 \v\-^ 
3 A2 



(9.9) 
(9.10) 



In Fig. 11 we have shown the configurations. 

In a similar way we can construct a more faceted height profile in terms of a gas of appropriately paired solitons 
in the slope field u with the only requirement that i) the solitons are well-separated so that they constitute saddle 
point solutions and ii) they satisfy periodic boundary conditions. For example a growing tip or the filling in of an 
indentation is described by the three-soliton configurations. A growing plateau formed by two step corresponds to a 
four-soliton configuration. We also notice that the two-soliton configurations corresponding to a moving step can be 
"renormalized" by the excitation of further two-soliton configurations corresponding to curvature of the step. In Fig. 
12 we have depicted the above special configurations. In Fig. 13 we have shown a general profile. 

X. "QUANTUM DESCRIPTION" OF A GROWING INTERFACE - FLUCTUATIONS 



In the previous section we demonstrated that the dominant morphology of a growing interface governed by the 
noisy Burgers equation in the weak noise limit can be described in terms of a dilute gas of propagating solitons. In the 
path integral the soliton contributions correspond to the stationary saddle poin ts in the {u, ip) phase space deter mine d 
by the principle of least action. By inspection of the one dimensional integral ( |6.l| ) and the saddle point result ( |6.2| ), 
it is clear that the soliton solution corresponds to the stationary point uq and the associated soliton action to S{uq). 
Consequently, we have not included the Gaussian fluctuations about the stationary point, yielding the multiplicative 
factor in Eq. ( |6.2|) of order A^/^ and depending on the second order derivative S"'(uo) evaluated at the stationary 
point, but only taken into account the exponential contribution determined by the action. In the context of the path 
integral the Gaussian fluctuations about the stationary soliton correspond to the linear diffusive mode spectrum in 
the presence of the soliton configurations and remains to be discussed. 

In order to proceed in the analysis of the path integral representation of the noisy Burgers equation we shall take a 
heuristic point of view and extract some information and physical insight by making use of the Feynman path integral 
structure of Z in Eqs. ( [3.2| - ^.4D ), deferring an analysis of the path integral per se to another context. The idea is to 
in a certain sense "deconstruct" the path integral and determine the form of the underlying "quantum field theory" 
leading to Z by the usual Feynman method 1 72 T^Jt^ . Since the slope field u and the noise field f in the path integral 
form a canonically conjugate pair with Poisson bracket (|]^), where u plays the role of a canonical "momentum" and 
ip a canonically conjugate "coordinate" , the first step, is to introduce the "quantum fields" u and ip satisfying the 
canonical commutator 



[u{x), (pix')] = —i — S{x — x') 



(10.1) 
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Here the ratio of the noise to the damping, A/z^, plays the role of an effective Planck constant just as in the path 
integral. We thus have an effective "correspondence principle" operating relating the "classical" Poisson bracket 
{A,B} to the "quantum" commutator [j4, J5], according to the prescription [A, -B] = —i{A/i'){A,B}. In a similar 
way the effective "quantum Hamiltonian" H is inferred from (5.8) [pO[, 



H = 



dx . 



(10.2) 



Whereas the fields u and (p by construction are Hermitian the Hamiltonian H is in general a non-Hermitian operator. 
Expressing H in the form Hq+Hi it is composed of an anti-Hermitian harmonic component Hq , governing the dynamics 
of the linear diffusive modes and a nonlinear Hermitian component iJi, describing the growth characterized by A. 
In the Heisenberg picture H is the generator of time translations and we obtain the usual Heisenberg equations of 
motion |81 



du .u ^ 



yielding "quantum field equations" of the same form as the "classical" field equations (6^6. £) 

dii 

'dt 

dip 

'at 



= —iiy'V ip + Xu'S/u 
= +wV^M + Xu\7ip . 



In a similar way, the momentum operator P, the generator of translation, is given by 



giving rise to the commutator relations 



V^ = z£[P,<p] • 



(10.3) 
(10.4) 

(10.5) 
(10.6) 

(10.7) 

(10.8) 
(10.9) 



Finally, the symmetry algebra in section V also holds in the "quantum case" by simply replacing the Poisson brackets 

by commutators according to the above "corr espondence principle" . 

The "quantum field equations" (10.5- 10^ )) together with the appropriate states of the Hamiltonian ( 10. 2| ) are 
completely equivalent to the path integral and thus provide an alternative description of the noisy Burgers equation. 
The noise-induced fluctuations in the slope field, represented by the different configurations or paths in the path 
integral weighted by the "classical" action 5, are replaced by "quantum fluctuations" in the underlying "quantum 
field theory" , resulting from the operator structure and the associated commutator algebra. We also note that the 
"quantum description" presented here is precisely the same as the one obtained in L based on the mapping of a 
solid-on-solid model to a continuum spin chain model in the quasi-classical limit. 



A. The Edwards- Wilkinson equation 



In order to demonstrate how the "quantum scheme" works it is instructive to evaluate the slope correlation function 
{uu){kuj) in Eq. (2.8) for the Edwards- Wilkinson equation. The dynamics of the EW case is governed by the 



unperturbed part of H in Eq. (10.2) 



Ho 



(V(^)^ 



dx 



(10.10) 



Introducing the usual "second quantization" scheme [p3[ in terms of Bose annihilation and creation operators at and 
aj. satisfying the commutator algebra [a^, aj^] — 6kk', we have for the slope and noise fields 
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2vL 



E 



ikx ] 1 



2vL 



yielding the unperturbed Hamiltonian and associated diffusive dispersion law 



-ivk 



(fO.ff) 
(f0.f2) 

(f0.f3) 
(f0.f4) 



Noting that the particle vacuum state |0) corresponds to a stationary state with average vanishing slope (0|m|0), i.e., 
a horizontal interface, it is an easy task to evaluate the slope correlation function. Since the path integral defines a 
Bose time-ordering [[73[ and using the time evolution operator, we have the identification 



{u{x,t)u{0,Q)) = (0|Tw(x,i)u(0,0)|0) = {Q\u{x)e-'"°\'\l^^/''^u{Q)\Q) . 
Using that evolves in time according to 

ak{t) = afc(0)exp ^-i^Lj"< 

we then obtain in Fourier space. 



2v 



{u{k, uj)u{—k, —Lo)) = i 

or in reduced form in complete agreement with Eq. (| 

{u{k, uj)u{—k, —uj)) 



|afc4|0) , (0|afc4|0)' 



Afc2 



(10.15) 



(10.16) 



(10.17) 



(10.18) 



This simple calculation demonstrates how the "quantum fluctuations" as expressed by the commutator algebra and 
the effective Planck constant A/i/ combine to produce the factor A in the correlation function which "classically" in 
terms of the EW Langevin equation originates from averaging over the noise Vry. 



B. The "quantum soliton" 



In the nonlinear case the "quantum dynamics" is governed by _ff = Hq + Hi in Eq. ( 10. 2| ). Introducing the Bose 
field ij} in configuration space. 



■0(x) — (1/VZ) Qfc exp [ikx) 



the Hamiltonian takes the form 

H = {~i){/^/v) / dx|V?/;t|2 - A(A/2i/)3/2 / ^^(^t „ ^)2v(v3t + ?/,) 



(10.19) 



(10.20) 



describing the many-body interaction between the linear diffusive modes governed by the first term Ho. 

Imposing the constraint of a horizontal interface we obtain (u) oc {ip — "0^) — 0, which implies that = (ip^)- 
Since the interaction term Hi does not conserve the number of particles, this constraint can only be satisfied for 
non-vanishing if the diffusive modes condense into a coherent condensate so that (ip) = {ip'') ^ 0. The resulting 
macroscopic wave function or condensate corresponds to the "classical" soliton mode discussed in the previous sections. 
The situation is quite similar to the phenomenological theory for supcrfluid Helium based on a condensate wave 
function. The condensate has two components, (^/i) and ('0^) or (u) ~ u and {(p) ~ and satisfies the coupled field 
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equations (6.5-6^), obtained from the "quantum field equations" (10.5-10.6) by ignoring "quantum fluctuatio ns" and 
replacing the terms, XuVu and AuVtp, by their average values, XuVu and XuVip. We can thus regard Eqs. (6.5-3.6) 
as two coupled Gross-Pitaevsky-type equations for the condensate wave function or soliton mode . 

The "classical" soliton is lo calized in space and carries energy and momentum, depending on the boundary conditions 
according to the expressions (B.2-8.3). Subject to "quantization" this mode becomes a bona fide "quantum mechanical" 

however, that the "quantum soliton" is delocalized 
A/v; here Aa;o is the uncertainty in the center of 

with well-defined 



quasi-particle with the same energy and momentum. Notice, 
owing to the "uncertainty principle" which implies that AxqAP ■ 

mass position for the soliton and AP the uncertainty in its momentum. For a "quantum soliton" 
momentum P and energy E we can in the usual way associate a wave number K and a frequency $7, according to the 
"de Broglie" relations, P = {A/i/)K and E = {A/i/)^}, and describe the quasi-particle by means of the wave function 
^' (X exp [—iQ,t + iKx]. Considering in particular a pair o f "quan tum solitons", describing a propagating step in the 
height profile with energy and momentum given by Eqs. (9.8-9.9), we obtain 



step 



/ v \ v\v\ 

= - (a) 4^ 



(10.21) 
(10.22) 



and the wave function takes the form 



^ oc exp [—iQt + iKx] = cxp [const(a:: 



Vph 



t)] 



(10.23) 



corresponding to a propagation with phase velocity Vpii 
for a localized soliton is the wave packet construction. 



(2/3)w. Noting, however, that the appropriate wave function 



"^WP oc y. exp [—iU,t + iKx] 



K 



obtained from a superposition of plane waves, we obtain the group velocity 



dK 



dQ/dv 
dK/dv 



(10.24) 



(10.25) 



This shows that the quasi-particle wave packet propagates with the same velocity as the "classical" soliton in complete 
accordance with "the correspondence principle" . Whereas the propagation velocity v determines the kinetics of the 
"classical" soliton, the energy and momentum are the fundamental characteristics in the "quantum" c ase; th e velocity 



V becoming the group velocity of the wave packet. We also notice from the wave packet form in Eq. ( 10.24) that the 



"quantum soliton" corresponds to a propagating mode. Finally, eliminating the velocity from Eqs. (10.21 
derive the "quantum soliton" dispersion law 



10.22) we 



n 



K 



where sign(Imii') — — signw. 



\K\ 



(10.26) 



C. "Quantum fluctuations" 



The final issue to consider in the qualitative "quantization" of the soliton system is the role of "quantum fiuctuations" 
in the presence of a "quantum soliton" . This problem is treated here by expanding the fields u an d i}> a bout a soliton 
or condensate configuration (uo, ^po). Inserting u — uq -\- 5u and (p = ipo + Stp in Eqs. (10.5) and (10.6) we obtain to 
linear order two coupled equations for the "quantum fluctuations" 6u and 5ip, 



dSu 

~dt 

d6(p 

"dt 



-iv\7'^5ip + AuoV^ii -t- A(Vmo)5<p 



(10.27) 
(10.28) 



These equations have the same form as the ones obtained by expanding in the Gaussian fluctuations about the 
stationary soliton solution in the path integral. 
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The equations of motion (10.27 jlb.28| ) describe the interaction of the linear "quantum diffusive modes" {Su,5(p) 
with the soUton configuration [uq, ipo) and constitute a generahzation to the noisy case of the linear stability equation 
in the analysis in A; the soliton again acts like a potential giving rise to phase shift effects and a gap in the diffusive 
spectrum. 



Like in the noiseless case, the equations (10.27-10.28) admit an analytical solution. Since the equations are Galilean 
invariant we need only consider the case of a static soliton. First noting that the soliton solution according to Eq. 
( [7.9[ ) is confined to the diagonal lines V(^o = iM^^^o, M — ±1, in phase space, the fluctuations Su, S(p are disentangled 
by transforming to "normal coordinates" along and perpendicular to the "soliton lines". Thus 1) introducing "normal 
coordinates" 6X ^ Sii + iSip, 5Y ^ Sit — i6ip, 2) inserting the static solution (7.11), uq = tanh (/csx), Vwo = 

/Li|u+|fcs cosh~^ (ksx), V(/9o = fJ-^uo, and 3) performing the scaling transformations SX h6X, 6Y h^^SX, where 
h = cosh^ (ksx), in order to absorb the linear terms in V, i.e., 



SX = 
6Y = 
h = 



h ^(Sii + iSip) 
h{5u — i5(p) 
cosh'^ {ksx) 

2v 



we arrive at the effectively decoupled equations for the "normal coordinates" 



d^>X 
dt 

dSY 



= +DSX + l)iyk^JY 
= -D5X + {^i+l)i'klSY 



(10.29) 
(10.30) 
(10.31) 

(10.32) 



(10.33) 
(10.34) 



Here the Schrodinger operator D has the same form as the stability matrix for the noiseless Burgers equation in A or 
the sine Gordon equation [l62 65|, 



cosh (kgx) 



(10.35) 



The wave number kg = X\u+\/ 2h' introduce d in section III depending on the soliton amplitude sets the inverse length 
scale. We also note that Eqs. (10.33 10.34) reduce to the linear case for A = since D —vV^ and /i — > 1. 

Since the Bargman potential cosh^ (ksx) admits an exact solution the spectrum of D defined by the eigenvalue 
equation D^n = i^^n^n is well-known [ p5[ and is discussed in A. It is composed of a zero-eigenvalue localized bound 
state mode and a band of phase shifted scattering modes. 



1 



" cosh(fcsx) 
Wo = 



^I/fc oc exp (ikx) 



k + iks tanh {ksx) 



ika 



Uk = -w{k^ + kl) . 
Expanding 5X and 5Y on a set of eigenstates 

SX = ^a„*, 
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(10.36) 
(10.37) 
(10.38) 
(10.39) 

(10.40) 
(10.41) 



we finally obtain equations of motion for the expansion coefficients a„ and 6„, 

IT 

dh 



at 



(10.42) 
(10.43) 



which we proceed to discuss. 
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1. The translation modes 



The zero- frequency of the Schrodinger operator D in Eq. ( |l0.35 ) is associated with the translation and bo osting 
of th e static sohton profile (uo,</Jo)- This is seen in the followi ng w ay. Since the "quantum field equations" (10.5- 
10.6 ) have the same form as the "classical" field equations ( |6!^|6.6[) they are equally satisfied by a soli ton solution . 
Co nsequently, a variation of the static soliton profile, ((5wo, is a solution of the linearized equations (10.27-10.28) 
or ( 10. 33| - 10.34 ) in the static case, corresponding to the bound state = 0. Furthermore, since the soliton depends 
parametrically on the center of mass position xq it follows that the fluctuations {5uq^5lpq) are proportional to the 
derivatives (Vwq, V(^o) with respect to xq, corresponding to a displacement of the soliton position. This mode is 
thus a translation or Goldstone mode associated with the broken translational symmetry and is a well-known feature 
of symmetry breaking "excitations" ; a Goldstone mode is excited in order to restore the broken symmetry |]3^ . A 
similar translation mode was also encountered in our discussion of the noiseless case in A. 

Focussing for example on the right hand static soliton for ^ 
the expansion coefficients 



1 and solving Eqs. (10.42-10.43) for n = we have 



ao = est. 



where 6q is the initial value for t — Q, and we obtain, using that Vuq oc cosh ^(fc^x), the fluctuation mode 

5u — ao + (6q -I- aoi'k'^t)V Uf) . 



(10.44) 
(10.45) 



(10.46) 



For ao = this mode corresponds to an inflnitesimal translation Sxq = &[] of the soliton, i.e., a change of the center 
of mass coordinate; for oq ^ the mode is equivalent to a boost of the soliton to a small velocity cx vkg. A similar 
discussion applies to Sip. 



2. The diffusive scattering modes 



The band of diffu sive scatteri ng modes are also easily discussed. From the equations of motion for the expansion 
coefficients in Eqs. ( 10.42 - 10.43 ) and again considering a right hand static soliton for fi — 1 we obtain ioi n ~ k the 
solution 



ak 



h = al 



-iujkt 



(10.47) 
(10.48) 



Here a° and 6^. are the initial values and the spectrum uj^ given by Eq. ( 10.39| ). For the fluctuation Sii we then have 

5u — ^^[flfc cosh (fcgx) -f hk cosh^"'^ {kgx)]^!, , (10.49) 



where ^'fc is given by Eq. ( 10.38|) . We note that 5u in the soliton case again is composed of both positive and 
negative frequency parts, exp {v{k^ + fc^)t) and exp {—v{k'^ -f kg)t), corresponding to growing and decaying modes in 
the stationary dr iven r egime, exhibiting a gap vk^ in the spectrum. In the linear EW case fcs = and Su assumes 
the form in Eq. (10.11). We shall not here dwell on the somewhat complicated ^-dependence but only observe that 
the main effect of the soliton on the diffusive modes apart from phase shift effects and spatial modulations is to lift 
the spectrum and create a gap vk^ = A^|u+p/4z^ depending on the soliton amplitude u+, the coupling A and the 
frequency v. 



D. Many-body description of a growing interface 



The above discussion of the "quantum solitons" and the "quantum diffusive modes" allow a heuristic qualitative 
discussion of a growing interface. The stochastic dynamics of the noisy Burgers equation (4.1) in the stationary regime 
can be rigorously interpreted in terms of a dilute Landau-type "quantum" quasi-particle gas composes of elementary 
excitations of two types: "Quantum solitons" and "quantum diffusive modes" . The "quantum mechanics" being 
equivalent to a Master equation description MjHI is basically relaxational corresponding to a complex Hamiltonian. 
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The elementary excitations fall in two classes: Linear diffusive modes and nonlinear soliton modes. 1) The linear 
diffusive modes are are associated with the damping term in the Burgers equation or, equivalently, the harmonic 
anti-Hermitian part in the Hamilto nian. T hese modes account for the relaxational aspects of the interface and are 
characterized by the dispersion law ( 10.39| ), i.e., 



Uk = -ii^ik"^ + k1) (10.50) 

with a gap vk^. As in our discussion of the noiseless case in A the gap can be associated with a non- vanishing current 
towards the center of the soliton where the damping in enhanced. 2) The nonlinear soliton modes are related to the 
nonlinear growth term in the Burgers equation or, equivalently, the nonlinear Hcrmitian part of the Hamiltonian. For 



a pair of solitons representing a growing step the dispersion law is given by Eq.(10.26), i.e., 



^K = ^-^[-) \K\-. (10.51) 

The mode is gapless and characterized by the fractional exponent 3/2. The soliton mode accounts for the growth 
aspects of the driven interface. For v ^ oo the linear damping dominates the growth and — > 0, also for A — + we 
attain the linear EW case and Six 0; finally for A ^ the stochastic aspects are quenched, solitons (and diffusive 
modes) are not kinetically or stochastically excited and U,k 0. 

In the linear EW case the fluctuating interface is in equilibrium and here described as a non-interacting gas of 
linear gapless diffusive modes. The statistical fluctuations appear as "quantum fluctuations" of the quasi-particle 
modes. Since the dispersion is quadratic we can also envisage the EW case as a "quantum" gas of free particles with 
imaginary mass. 

In the nonlinear Burgers case the "quantum soliton" emerges as a new additional quasi-particle, corresponding to 
the faceted genuine growth of an interface. The linear modes become subdominant in the sense that they develop a 
gap in the spectrum and correspond to superposed damped "ripple modes" on the soliton configurations. The diffusive 
modes extend over the whole configuration and are phase-shifted due to refiectionless scattering off the solitons like 
in the noiseless case. A scattering analysis also, in accordance with Levinson's theorem, shows the diffusive spectrum 
is depleted by a number of states, corresponding to the translation modes of the solitons. 

Before turning to the heuristic scaling analysis in the next section, we wish to add a few more remarks concerning 
the structure of a field theoretic or many-body description of the noisy Burgers equation. There are basically two 
equivalent modes of approach: 1) A direct evaluation of the path integral in the weak noise limit in a saddle point 
approximation for a dilute soliton gas, including the diffusive modes, corresponding to Gaussian fluctuations about 
the saddle points, and summing over periodic orbits in order to include secular effects or 2) A construction of the 
equivalent "quantum many-body theory" on the basis of the "quantum representation" of the path integral. Including 
the soliton modes as space- and time-dependent condensate configurationsi, as mentioned in section XB, the many- 
body approach is similar to the microscopic theory on interacting bosons ||8^ , |9l| , |9^ with anomalous propagators etc. 
There are, however, some notable differences. In the case of interacting bosons the uniform condensate acts as a 
particle reservoir and changes the free boson dispersion law a; oc to a linear acoustic phonon branch uj oc p. In 
the present case , the condensate corresponding to a soliton or gas of solitons is non- uniform and time-dependent, 
governed by the "classical" field equations. The free diffusive modes with dispersion cx develop a gap depending 
on the soliton amplitude or velocity, whereas the soliton or condensate mode emerges as a new quasi-particle with 
dispersion u) cx p^^'^. 

XI. SCALING AND UNIVERSALITY CLASSES 

In addition to providing a many-body description of a growing interface in terms of a Landau-type quantum quasi- 
particle gas of propagating "quantum solitons" and damped "quantum diffusive modes" , the path integral formulation 
also offers as a by-product some insight into the scaling properties, i.e., the behavior of the interface in the limit of 
large distances and long times. 

We shall here focus on the scaling properties of the slope correlation function summarized in the dynamical scaling 



form (1.7), i.e., assuming t = in the stationary regime, 

(u(x,i)"(o,o))-ixp(«-i)/(|i|/|xr) (11.1) 

The scaling issue is then to determine the roughness or wandering exponent C, the dynamic exponent z, and the 
scaling function f{w). 
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In the EW case the scaUng function / is given by (2.9), i.e., 



f{w) = {A/2iy){4:TTiy)-^/^ur'^/'^exp[-l/4:i'w] , 



(11.2) 



implying the exponents (C, z) — (1/2,2). In the Burgers ca se Ga Ulean invariance leads to the scaling law (1.10), i.e., 
^ + z = 2, which together with the stationary distribution ( 1.11 ), an effective fluctuation-dissipation theorem, yields 
the exponents {C,z) = (1/2,3/2). 

According to the path integral formulation in sections IV and V, using Eqs. (4/7) and (|5.2|) the slope correlation 
function is given by 



(M(a;,t)M(0,0)) = Z{Oy 



dudip 



exp 



. A . 



u{xt)u(00) 



(11.3) 



or in terms of the underlying "quantum field theory" , noting that the path integral by construction defines time- 
ordered products fz^, 



{u{x,t)u{0,0)) = (0|T{t(a;,i)u(0,0)|0) 



(11.4) 



Here |0) denotes the appropriate stationary st ate f or the system |86[ . 

In order to elucidate the structure of Eq. ( |11.4[ ) we construct a spectral representation by i) displacing the slope 
field u{x, t) to the origin in space a nd tim e by means of the Hamiltonian H and the momentum P, using the integrated 
form of the commutator relations ( |l0.3|) and ( |10.8[ ), and ii) inserting intermediate eigenstates \P) with momentum P 
and energy Ep. The first step implies the relation 



u{x, t) = exp 



i—{Px + Ht)^u{0,0) exp 



-i-^{Px + Ht) 



(11.5) 



secondly, inserting intermediate states, using H\P) = Ep\P) and P\P) — P\P): introducing the frequency and wave 
number {^Ik^K) associated with the elementary excitations or quasi-particles, and lumping the matrix elements in 
an effective form factor, G{K) ~ {{)\u\K) {K\u\Q) , we arrive at the spectral representation 



dK 

{u{x,t)u{Q,Q)) ^ / —G{K)eyiY,[-i{^K\t\-iK\x 
' 27r 



(11.6) 

X, imply evenness in the dependence on x and 



The time-ordering in Eq. (11.4) together with parity invariance, x - 
t; also G{K) must be even in K . 

The spectral form (11.6) is only schematic. For a multi-soliton diffusive mode intermediate eigenstate {fcj}), 
where Ki and kj denote the soliton and diffusive mode wave numbers, respectively, with total wave number K = 
Ki + kj and total frequency il — Q.Ki + , we have strictly speaking the spectral form, say for t > 0, 



{u{x,t)u{0,0)) 



J []dif,dfc^.G({i^,},{%})e-*'^^-'''+^.'^^"+^^-''"'+^."'=^^*l . (11.7) 



Since the soliton are transparent with respect to the diffusive modes as discussed in section XC, the operator ii only 
excites a single mode k extending across the system, i.e., G{{Ki}, {kj}) ~ G{{Ki}, k), and assuming furthermore that 
G{{Ki], k) factorizes approximately in accordance with the dilute soliton gas picture, G{{Ki\^ k) ^ Goik) Hi Gs{Ki), 
we obtain summing over the solitons 



(w(x,t)w(0,0)) 



JdkGoikje 



-ikx — iujt 



\~ j dKGs{K)e- 



(11., 



The expression (11.8) is clearly not correct in detail since we have not solved the many-body problem but only made 
some plausible assump tions concerning the form factor G. Nevertheless, from the point of view of discussing the 
s calin g properties Eq. ( 11. S| ) has th e req uired structure and serves our purpose. In the EW case Gs{K) = and Eq. 
( 11.8 ) reduces to the scaling form ( 11. 2| ). In the Burgers case J7 oc ilXp^, i.e., exp (— ifit) = exp (-hconst.|X|3/2i), 
and the denominator (1 — / dKG s{K)e^^^^^''^^^)^^ ^ f d KG s{K)e~^^^~^^*' controls the scaling behavior. In both 
cases we can use the simplified general spectral form (11.6). 

For the purpose of a discussion of the scaling properties we first consider a general quasi-particle dispersion law 
with a gap A, stiffness constant A, and exponent /3, 
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nK = A + A\Kf 



(11.9) 



For large distances \x\ ^ a, where a is a microscopic length defining the UV cut-off K ^ 1/a implied in Eq. ( 11.6 ), the 
spectral representation samples the small wave number regions K <^ 1/a. Assuming that the form factor is regular for 
s mall K, i.e., G{K) ~ G(0) + (l/2)i4r^G"(0), and rescaling Kx — > we obtain, inserting the general dispersion 
(11.9), the spectral representation in a more appropriate scaling form 



(7.(x,t)7.(0,0))~G(0)e-'^*x-i / !iiie-^^l^l"l*l/ 



dK 



iK 



2-K 



(11.10) 



We emphasize again that the spectral representation (11.10) can only be considered as a heuristic expression since 
we have not here carried out a detailed analysis of the non-Hermitian non-Lagrangian field theory underlying the path 
integral. Nevertheless, we believe that we can already here draw some interesting general consequences concerning 
the scaling properties of a growing interface. 

First we observe that in the presence of a gap A 7^ there is no scaling behavior. For the diffusive mode in the 
presence of a soliton the spectrum is for example given by Eq. (10.50) with a gap A — —ivk'^, implying an exponential 
fall-off with t in Eq. ( |ll.lC| ). Consequently, only gapless excitations for A = contribute to the scaling behavior. 
The gapless excitations are associated with the so-called zero temperature fixed point behavior of the "quantum field 
theory" and determine the scaling properties. 

11.10|) for A = with the scaling form ( 11.1 ) we immediately identify 
1/2 and the dynamic exponent z — (3. We also note that whereas the exponent C = 1/2 



Furthermore, comparing the spectral form 
the roughness exponent C, 



essentially follows from a simple regularity property of the form factor G{K) with leading term G(0) for small the 
exponent z is tied to the exponent (3 in the quasi-particle dispersion law. 

In the linear EW case the diffusive gapless modes with dispersion law (10.14), i.e., ujk = —ivk"^, exhaust the 
spectrum and we obtain (3 ^ z = 2, corresponding to the EW universality class in Table 1. Also the spectral form 
yields the scaling function ( 11.2 ) with the identification G(0) = A/2j/. 

In the nonlinear Burgers-KPZ case the soliton modes with gapless dispersion (10.51), i.e., oc {/S./vY/'^\\K\^/'^ ^ 
exhaust the bottom of the spectrum and yields the exponent (3 = z — 3/2, corres pondin g to the Burgers-KPZ 
universality class in Table 1; the linear diffusive modes develop a gap according to Eq. ( 10.50| ), become subdominant 
and do not contribute to the scaling behavior. 

The above discussion thus provides a dynamical interpretation of the scaling properties, exponents, and universality 
classes. The universality class is determined by the lowest-lying gapless excitation. The spectral form also elucidates 
the robustness of the roughness exponent C which is the same for both universality classes. In the case of the stationary 
equal-time fluctuations we set t = Q va Eq. ( 11. §| ) and the resulting scaling form yielding C, does not depend on the 
specific quasi-particle dispersion law; this argument is equivalent to the effective fluctuation-dissipation theorem 



(11. 



yielding the stationar y distr ibution (1.11) independent of the nonlinear drive A 
spectra] 
with Eq 



The spectral form (11. IC) also provides an expression for the scaling function f{w). 
~ (PI) 



f{w) = G(0) 



dK 



i(A\K\'w+K) 



We obtain comparing Eq. 



(11.11) 



The above expression is at best heuristic but we do notice that it has the correct limiting behavior, i.e., f{w) ~ const 
for w — > and f{w) ^ w^^^^ for w 00. 

The scaling function f{w) describing the behavior of the strong coupling fixed point has been accessed both 
numerically [|87|-p9| and by means of an analytical mode coupling approach [EslpQl , based on a self-consistent one-loop 
calculation, i.e., to first order in A, and assuming vanishing vertex corrections. The agreement between the numerical 
simulations and the analytical method is good, indicating that the mode coupling approach seems to capture essential 
properties of the strong coupling fixed point behavior. 

The heuristic and preliminary character of the scaling function given here does not allow a detailed comparison. We 
note, however, that since A cx A(A/i/)^/^ in the Burgers-KPZ case the dimensionless argument in the scaling function 
f{w) is A((A/i/)^/^t/a:'^/^. This is in complete agreement with the driven lattice gas DRG analysis in and with 
the general arguments advanced in the mode coupling analysis in |^ , ^ . 

We also note the curious fact that the spectral form (11.6) for a gapless dispersion with exponent (3 bears resemblance 
to the form of the probability distribution for a one-dimensional Levy fiight with index = (3 . The case fi = f3 = 2 
corresponds to ordinary Brownian walk, whereas /x = /? = 3/2<2is equiva lent to super diffusion. 

The present analysis of the scaling properties based on the spectral form (11.6) originates from a weak noise saddle 
point approximation to the path integral and as such only holds for A ^ 0. However, within the general assumptions 
underlying the application of scaling theory and the notion of universality classes, we expect the exponents and scaling 
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function to be universal characteristics of the system and thus independent of the noise strength A. This property can, 
however, be reconciled with the present many-body approach if we assume that an enhancement of the noise strength, 
that is a stronger drive of the system, only leads to a dressing of the quasi-particle spectrum, i.e., a change in the 
stiffness constant, and not to a change in the exponent /3. In the "quantum mechanical" language this corresponds 
to the assumption that the WKB approximation also holds in the strong "quantum regime" as far as the exponent of 
the quasi-particle dispersion law is concerned. 

We conclude this section with a few speculative remarks concerning the "breakdown of hydrodynamics" . The noisy 
Burgers equation is basically a nonlinear conserved hydrodynamical equation derived by combining the conservation 
law du/dt + \7u = with a constitutive equation for the current, j = ~v\7u — {\/2)v? — ry, with transport coefficients 
v (the damping) and A (the mode coupling). The expression for the deterministic part of j is thus based on a gradient 
expansion to lowest order and the simplest quadratic nonlincarity in u. The issue is in which way the mode coupling 
term affects the hydrodynamical properties. In this context "breakdown of hydrodynamics" usually refers to the 
situation where the underlying regularity structure of the gradient expansion, i.e., in wave number space regularity 
in an expansion in fc, breaks down. 



In the present many-body formulation, entailing the spectral form ( 11.6 ), we obtain in frequency space 



(uu)(fc,u;) - Re (11-12) 

LO — \lk 



In the linear EW case for uj^ — —ivk'^ we recover the diffusive form (2^), corresponding to a diffusive pole oj^ ~ —ivk^ 
in the complex lu - plane. However, in the nonlinear mode coupling case for A 7^ 0, fifc cx |fcp/^ and [uu) (fc, to) develops 
a branch cut structure, corresponding to a non-analytic wave number dependence in the current, i.e., a breakdown of 
hydrodynamics. 

XII. DISCUSSION AND CONCLUSION 

In the present paper we have advanced a novel approach to the growth morphology and scaling behavior of the 
noisy Burgers equation in one dimension. Using the Martin-Siggia-Rose (MSR) technique in a canonical form we 
have demonstrated that the physics of the so far elusive strong coupling fixed point is associated with an essential 
singularity in the noise strength and can be accessed by appropriate theoretical soliton techniques. 

The canonical representation of the MSR functional integral in terms of a Feynman phase space path integral with 
a complex Hamiltonian identifies the noise strength as the relevant small non-perturbative parameter and allows for 
a principle of least action. In the asymptotic weak noise limit the leading contributions to the path integral are given 
by a dilute gas of solitons with superposed linear diffusive modes. The canonical variables are the local slope of the 
interface and an associated "conjugate" noise field, characteristic of the MSR formalism. In terms of the local slope 
the soliton and diffusive mode picture provide a many-body description of a growing interface governed by the noisy 
Burgers equation. The noise-induced slope fluctuations are here represented by the various paths or configurations 
contributing to the path integral. 

The canonical formulation of the path integral and the associated principle of least action also allow us to associate 
energy, momentum, and action with a given soliton configuration or growth morphology. This gives rise to a dynamical 
selection criterion similar to the role of the Boltzmann factor exp [—E/T] in equilibrium statistical mechanics which 
associates an energy E with a given configuration contributing to the partition function; in the dynamical case 
the action S provides the weight function for the dynamical configuration. More detailed, in the dynamical case, 
"rotating" the noise variable, ip — > — «<^, we have the partition function 

ZdynC^ / ]^(iMd(pexp [--^^J (12.1) 

dxdt[u—-H{u,ip)] (12.2) 

n^-'^[{\7uf-{\7^f] + ^u'\/p, (12.3) 
whereas in the equilibrium case we have the general form in one dimension, 

Zeq(x lY[dpdqexp[-^H] (12.4) 

X 

H = j dxn{p,q) , (12.5) 
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where H is the Hamiltonian density. By comparison we note that the noise strength A plays the role of a "temperature" 
in the dynamical case. We also observe that Zdyn for the dynamical 1-D problem, treating time as an additional 
coordinate, i.e., t ^ y, is equivalent to a 2-D equilibrium partition function with Hamiltonian 



and temperature v / /S.. 

In addition to providing a physical many-body picture of the morphology of a growing interface in terms of soliton 
modes accounting for the growth aspects and diffusive modes corresponding to the relaxational aspects, the present 
approach also gives insight into the scaling properties. The perspective here is not a "coarse-graining" procedure, 
replacing the original description by a scaling description with ensuing dynamical renormalization group (DRG) 
equations, but rather a focus on the gapless elementary excitations or quasi-particlcs of the many-body theory. 

The case of simple scaling characterized by a roughness exponent, a dynamical exponent, and a scaling function, 
corresponding to a simple fixed point structure in the DRG analysis, is here represented by a simple quasi-particle 
mode exhausting the bottom of the spectrum with a gapless dispersion law. The dynamic exponent is given by the 
exponent in the quasi-particle dispersion law, whereas the roughness exponent follows from a regularity property of 
the form factor in a spectral representation of the slope correlations; the scaling function being given by the spectral 
form itself. 

Our analysis shows that the nonequilibrium growth dynamics in one dimension is controlled by solitons or dynamic 
domain wall. In this respect their is a parallel between the present kinetic growth problem and other well-studied 
low dimensional equilibrium problems also controlled by localized excitations such as the one dimensional Ising model 
with domain wall excitations or the two dimensional XY model characterized by vortex excitations ]36| . The present 
approach is conducted in one dimension and assumes a spatially short-range correlated conserved noise in order to 
implement the shift transformation leading to the canonical formulation and the separation of the Hamiltonian in a 
harmonic part and an interacting part. In higher dimension the Burgers equation becomes a vector equation with a 
nonlinear term X{u\J)u and we obtain a more complicated Hamiltonian governing the dynamics. We have not pursued 
the higher dimensional case yet, but it is nevertheless interesting to speculate that the strong coupling fixed point 
behavior in d = 2 and below is associated with higher dimensional localized soliton- like excitations, i.e., the dynamics 
is defect-dominated. 

Finally we wish to comment on some recent work on the driven Burgers equation with noise at large length scales 
modelling forced turbulence. This problem has been treated using a variety of methods such as operator product 
expansions , instanton calculations pst , and replica methods [ |9q | . In this context the non-perturbative instanton 
methods used in order to determine the tail of the velocity probability distribution seem superficially to be related 
to the present soliton approach in that they are also based on a saddle point approximation to the MSR functional 
integral. However, unlike the the soliton which is only localized in space, the instanton is localized in both space and 
time, yielding a finite action, and does not represent a soliton in the Burgers equation. The difference between the 
two approaches is related to the spatial correlations of the noise driving the system: In the forced turbulence case the 
noise is correlated at large distances, whereas in the growth problem the noise is assumed to be conserved and delta 
function correlated in space. 



Discussions with J. Krug, M. Kosterlitz, M. H. Jensen, T. Bohr, M. Howard, K. B. Lauritsen and A. Svane are 
gratefully acknowledged. 
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FIGURES 



u = Vh 



FIG. 1. We depict the general growth morphology for a 1-D interface in terms of the slope field u and the height field h 
(arbitrary units). 
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(a) 



|<uu> (k,(o) 






FIG. 2. In a) wc depict the slope correlation function for the diffusive mode in the linear EW case. The Lorcntzian is 
centered about lu = and has the "hydrodynamical" line width i^k^ vanishing in the long wave length limit. In b) we show the 
scaling functions / for the space and time-dependent slope correlations. For large w f falls off as w~^^^ , for small ui / — > 
with an essential singularity for w = 0. / is peaked at the value ~ 0.12A/i^ for w = l/2v. In b) we have shown the scaling 
function g for the wave number- frequency dependent correlations. / has a Lorentzian form with height A/i/ and width ~ v^^'^ 
(arbitrary units). 
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FIG. 3. We show a single moving soliton profile propagating to the left and the corresponding smoothed cusp in the growth 
profile. This configuration is driven by currents at the boundaries, corresponding to non-vanishing u± and is persistent in time 
(arbitrary units). 
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FIG. 4. We show the diffusive dispersion law in the presence of a soliton. The gap in the spectrum is given by uk, = 
where w+ is the soUton amphtude. The dashed hne indicates the gapless spectrum in the hnear case (arbitrary units). 
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FIG. 5. We here depict the transient evolution of the slope field u from an initial configuration uo in the case of the noiseless 
Burgers equation. We have also shown the evolution of the associated height field h. The transient morphology consists of 
propagating right hand solitons connected by ramps (arbitrary units). 
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FIG. 6. Here we depict in graphic form the basic principle of an asymptotic steepest descent or saddle point calculation. 
The leading contribution to the integral /(A) in Eq. (^^) arises from the saddle point uq (in the figure a minimum) and nearby 
fluctuations Su = uq — u (arbitrary units). 
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FIG. 7. In a) we show the "classical" path corresponding to the stationary point of the action S in the weak noise limit 
and nearby paths corresponding to fluctuations. In b) we show the saddle point regions in (u, if) phase space (arbitrary units). 
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FIG. 8. In a) and b) we show the static right and left hand solitons and the smoothed static downward and upward 
the associated height field (arbitrary units). 
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FIG. 9. We show the overlap between two well-separated solitons (arbitrary units). 
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FIG. 10. In a) and b) we show right and left-moving soUtons with vanishing amphtude at x = —L/2. The associated height 
profiles correspond to the bottom and top part of a step, respectively (arbitrary units).. 
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FIG. 11. In a) and b) we show two two-soliton configurations moving with opposite velocities corresponding to the left and 
right propagation of a step in the height profile (arbitrary units) . 
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FIG. 14. We depict the quadratic diffusive dispersion law with gap A^w+/4i/ and the gapless soUton dispersion law with 
fractional exponent 3/2 (arbitrary units). 
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